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ABSTRACT 

We use a set of simulation-based models for the dissipationless evolution of galaxies since z = 1 to constrain 
the fate of accreted satellites embedded in dark matter subhalos. These models assign stellar mass to dark 
matter halos at z = 1 by relating the observed galaxy stellar mass function (GSMF) to the halo+subhalo mass 
function monotonically. The evolution of the stellar mass content is then followed using halo merger trees 
extracted from iV-body simulations. Our models are differentiated only in the fate assigned to satellite galaxies 
once subhalos, within which satellites are embedded, disrupt. These models are confronted with the observed 
evolution in the massive end of the GSMF, the z ~ brightest cluster galaxy (BCG)-cluster mass relation, and 
the combined BCG and intracluster light (ICL) luminosity distribution — all observables expected to evolve 
approximately dissipationlessly since z = 1 . The combined observational constraints favor a model in which the 
vast majority (> 80%) of satellite stars from disrupted subhalos go into the ICL (operationally defined here as 
light below a surface brightness cut of /x, 23 mag arcsec" 2 ). Conversely, models that leave behind a significant 
population of satellite galaxies once the subhalo has disrupted are strongly disfavored, as are models that put a 
significant fraction of satellite stars into the BCG. Our results show that observations of the ICL provide useful 
and unique constraints on models of galaxy merging and the dissipationless evolution of galaxies in groups and 
clusters. 

Subject headings: galaxies: evolution — galaxies: halos — galaxies: mass function — galaxies: clusters — 
cosmology: theory — dark matter 



1. INTRODUCTION 

The formation and evolution of massive (M stal > 10 h Mq) 
elliptical galaxies is thought to be inexorably linked to the for- 
mation and evolution of the large scale structure of the Uni- 
verse. The classical picture where in massive elliptical galax- 
ies form "monolithically" at z > 5 (Partridge & Peebles ?! 9671) 
has been replaced by more nuanced scenarios that decouple 
the epoch at which these galaxies formed the bulk of their 
stars from the epoch (or epochs) at which these stars were 
assembled to form the final galaxy. These more complex sce- 
narios arise fairly naturally within the context of the hierarchi- 
cal grow th of structure in the now favored ACDM cosmology 
(see e.g. lBaugh et alJll996UNeistein et al.ll2006l) . While stel- 
lar population modeling has firmly place d the epoch of star 
formation in these galaxies at z > 2 (e.g. | Bower et alJI 1 992k 
Trager et al. 2000; van Dokkum & Fran xbOOU lThomas et al. 
20051: iJimenez e t al. 2006), the assembly history of massive 



galaxies is still far from clear, and is the focus of this work. 

The evolving space density of massive galaxies over time 
provides important clues to their assembly history. Evolution 
in the galaxy stellar mass function (GSMF) si nce z ~ 1 ap- 
pears quite mild for the most massive galaxies dFontana et all 
2004tlDrorv et alj2004trBundv et all2005UBorch et al Il2006h 
Fontana et a l. 2006; Cim atti et alj2006tlAndreonl2006l) . Esti- 



mates of evolution in the luminosity function are also consis- 
tent with massive galaxies pass i vely evolving from z ~ 1 to the 
prese n t (|Cirasuolo et alj|2005 IWake et al J 12 006: Fab er et all 
I2006t IWillmer et alJ l2006t iBrown et alJ 120071: ICaputi et alJ 
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2006; iBlantonl 120061) . These same observations also find 
roughly a doubling in the total stellar mass density from 
z^ltoz^O — the implication being that star formation 
occurs primarily in less massive galaxies at z < 1. In addi- 
tion to number counts, estimates of the merger rate of mas- 
sive galaxies can in principle constrain their assembly history, 
though current observations fail to present a consistent p icture 
jvan Dokkumll2005l; iBell et alj|2006t iMasiedi et ai]|2006l) . 

The majority of massive elliptical galaxies are (or have been 
in their recent past) the brightest cluster galaxies within large 
group- or cluster-sized halos (referred to simply as "clusters" 
in the remainder of the paper), located near the centers of the 
halo potential well. It is thus interesting to study the formation 
of such galaxies in the general context of cluster formation. If, 
as recent observations suggest, the majority of massive galax- 
ies were already in place at z ~ 1, then on the surface it ap- 
pears difficult to reconcile this with the much more substantial 
evolution of their host dark matter halos (massive halos grow 
by factors of > 3 in mass since z = 1). It is one of the goals of 
the present work to address and resolve this tension. 

Within the ACDM framework, groups and clusters of galax- 
ies are expected to be continually accreting new galaxies. Af- 
ter entering a cluster, the stars in satellite galaxies can 1) be all 
deposited onto the central galaxy, 2) stay bound as a satellite 
galaxy, or 3) be scattered into the intra-cluster light (ICL). For 
the purposes of this study we define the ICL as the stars be- 
yond the optical radius of the central galaxy, i.e. the light not 
accounted for in the photometry of the central galaxy itself. In 
reality, a combination of these possibilities can occur. For ex- 
ample, a satellite galaxy can be partially stripped and deposit 
a fraction of its stellar mass into the ICL before merging and 
contributing the rest of the stars to the central galaxy. In addi- 
tion, when two massive galaxies merge, a certain fraction of 
stars will acquire large kinetic energy and will move to large 
radii, outside the optical radius of the remnant. While evo- 
lution in the space density of massive galaxies strongly con- 
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strains the importance of the first scenario, this observation 
cannot readily distinguish between scenarios two and three. 

Fortunately, fates 2) and 3) have effects on observable prop- 
erties of the ICL, which can thus constrain the amount of stars 
that can be lost to the ICL during mergers or tidal stripping. 
These observations suggest that the fract ion of total cluster 
light bound up in the ICL is ~ 10-30% (IZibetti et al.ll2005l; 
IGonzalez et"aLll2005b iKricket al l [20061) . with the ICL com- 
prising a significant fraction (~ 50-8 0%) of the combined 
light from the cen tral galaxy and ICL (IGonzalez et alj|2005t 
ISeigar et alll2006L note that this fraction depends sensitively 
on the way that the ICL and BCG are separated). Addition- 
ally, models in which the stellar component of satellite halos 
is never disrupted by tides would have quite different predic- 
tions for the number of satellite galaxies in a halo of a given 
ma ss, and thus for the small-s cale clustering of galaxies (see 
e.g. lBerlind & Weinberg 2002). Indeed, the clustering of mas- 
sive galaxies in combination with their evolving spa ce den- 
sity has recently been exploited by Whi te et ail (12007) to con- 
strain the disruption rate of massive satellite galaxies between 
z~0.9andz~0.5. 

The goal of this study is to confront these and other ob- 
servational constraints simultaneously with a series of sim- 
ple, simulation-based models in order to gain insight into 
the fate(s) of satellite stars. The models presented herein 
combine a simple prescription for relating galaxies to dark 
matter halos at z ~ 1 with the assembly history of these 
halos extracted from A^-body simulations in order to fol- 
low the dissipationless growth of massive galaxies to z ~ 
0. The relation between galaxies and halos at z ~ 1 is 
generated by assigning the most massive galaxies to the 
most massive halos monotonically. This has been shown 
to successfully reproduce a wide variety of observations 
dColin et alj|1999t iKravtsov & Klypiril 1 1 9991: iNevrinck et all 
2004 iKravtsoy et alj 12004 IVale & Ostrikeri 12004 120061: 
Tasitsio mi et al. 20 04 IConrov et aTl 120061 : IShankar et alJ 
20061: IVale & Ostrikeri l2007h . Note that this model consid- 
ers both subhalos, which are halos contained within the virial 
radii larger halos, and what we will call distinct halos, which 
are halos not contained within the virial radii of larger halos. 
We follow the dynamical evolution of subhalos after they ac- 
crete onto their host halo using merger trees extracted directly 
from cosmological simulations, rather than a semi-analytic 
model. 

Semi-analytic models (SAMs) for the formation and evo- 
lution of galaxies within a cosmological context, depend- 
ing on the adopted assumptions about galax y formation 
physi c s, are capable of pred icting both strong (jBaugh et all 
1996; [De Lucia etal] 120061) and mild jBower et al] l2006t 
Kitzbichler & White! 120071: IMonaco et all 120061) evolution in 
the number density of massive galaxies since z ~ 1 . The most 
massive galaxies in many of these models have formed the 
bulk of their stars at z > 2, in agreement with observations. 
Hence differences between these models are due primarily to 
different treatments of the assembly history of the massive 
galaxies. 

Our approach is sim ilar in spirit to that of a recent study by 
IMonaco et alJ (|2006) who used a SAM to follow the evolu- 
tion of galaxies. These authors artificially turned off star for- 
mation at z < 1 in order to follow the dissipationless growth 
of galaxies at late times, similar to what we do here. The or- 
bital evolution of satellites in their SAM was computed with 
simple analytical approxim ations to dynamical f riction, tidal 
heading and tidal stripping. IMonaco et alJ ((2006) showed that 



the observed evolution in the space density of massive galax- 
ies is reproduced in their model only if they allow for > 30% 
of stars from disrupted satellites to be transfere d into the ICL. 
The pr esent work goes further than the study of lMonaco et alJ 
(2006) by 1) using an independent set of simulations with 
satellite tracks extracted directly from the simulations and 2) 
comparing to a wider array of observations and hence provid- 
ing more general constraints on the fates of the stars within 
satellite galaxies. 

The rest of this article unfolds as follows. In f|2]we describe 
the simulations, halos catalogs and merger trees used in this 
analysis. Sj3]ouflines the details of our models and ^contains 
comparisons between the models and several observations. 
The implications of these results and comparison to related 
work is discussed in Sj5] Throughout this paper we assume 
a ACDM cosmology with (tt m ,Q, A ,h, cr 8 ) = (0.3,0.7,0.7,0.9), 
except in ||2]where we leave quantities in terms of the reduced 
Hubble constant, h. 

2. SIMULATIONS, HALO CATALOGS, AND MERGER TREES 

The simulations used here were ru n with the Adaptive 
Refinement Tr ee (ART) JV-body code (IKravtsov etalj|1997t 
lKravtsovll99 9). The ART code implements successive refine- 
ments in both the spatial grid and temporal step in high den- 
sity environments. These simulations were run in the concor- 
dance flat ACDM cosmology with Q m = 0.3 = l-Q\, h = 0.7, 
where f2 m and are the present-day matter and vacuum den- 
sities in units of the critical density, and h is the Hubble pa- 
rameter in units of 100 km s" 1 Mpc" 1 . The power spectra used 
to generate the initial conditions for the simulations were de- 
termined from a direct Boltzmann code calculation (courtesy 
of Wayne Hu). We use a power spectrum normalization of 
as = 0.90, where erg is the rms fluctuation in spheres of 8/T 1 
Mpc comoving radius. 

The simulation used herein was run in a box of length 120 
h Mpc with particle mass m p = 1.07 x 10 9 h~ x M Q , peak 
force resolution of h pea k = 1 .8 h kpc, and 5 12 3 particles. We 
have checked that our results remain unchanged when utiliz- 
ing a smaller box with smaller particle mass (a box of length 
80 hr 1 Mpc with particle mass m p = 3.16 x 10 8 h~ l M Q ). 

From this simulation we generate dark matter halo catalogs 
and dark matter halo merger trees. Our models rely not only 
on distinct halos, i.e., halos with centers that do not lie within 
any larger virialized system, but also subhalos, which are lo- 
cated with the virial radii of larger systems. When we refer to 
a "halo" generically we mean both distinct halos and subha- 
los. 

Distinct halos and subhalos are identified using a variant 
of the bound densit y maxima (BDM) halo finding algorithm 
iKlvpin et al.1 d 19991) . Details of the a l gorith m and parameters 
used can be found in IKravtsov et al.1 j2004l) : we briefly sum- 
marize the main steps here. All particles are assigned a den- 
sity using the smooth algorithm 4 which uses a symmetric 
SPH smoothing kernel on the 32 nearest neighbors. Starting 
with the highest overdensity particle, we surround each poten- 
tial center by a sphere of radius rg n( j = 50h~ l kpc and exclude 
all particles within this sphere from further search. Hence no 
two halos can be separated by less than rg n( j. We then con- 
struct density, circular velocity, and velocity dispersion pro- 
files around each center, iteratively removing unbound par- 

4 To calculate the density we use the publicly available code smooth : 
http : / /www-hpec . astro. washingt on. edu /tools/ 
tools . html 
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tides as described in iKlvpin et al.l d 1999t) . Once unbound 
particles have been removed, we measure quantities such as 
V max = yjGM(< r)/r| max , the maximum circular velocity of 
the halo. For each distinct halo we calculate the virial radius, 
defined as the radius enclosing overdensity of 334 with re- 
spect to the mean density of the Universe at the epoch of the 
output. We use this virial radius to classify objects into dis- 
tinct halos and subhalos. The halo catalogs are complete for 
halos with more than 50 particles, which corresponds, for the 
box with length 120 h~ l Mpc, to 5.35 x lOVMg. 

Halo merger trees have also been constructed for this simu- 
lation (f or a detailed de scription of the merger tree construc- 
tion, see Allgood 2005). These merger trees allow us to tabu- 
late M^P? for subhalos, the virial mass at the time when a sub- 
halo first crosses the virial radius of a distinct halo. Since sub- 
halos are subject to dynamical processes such as tidal strip- 
ping, M^ r c will always be greater than or equal to the present 
M v j r . This accretion epoch quantity is used in our models, to 
which we now turn. 

3. THE MODELS 

3.1. Connecting galaxies to halos 

Recent studies have shown that models in which galax- 
ies are associated with the centers of dark matter halos and 
subhalos accurately reproduce a wide variety of observa- 
tions both at low and high redshift (e ■g.lKravtsov et al.ll2004t 
I Vale & Ostrikerl 12004 IConrov et all 120061) . Herein, stellar 
mass is assigned to each halo in a simulation by assuming a 
monotonic relation between stellar mass and halo virial mass 
using the observed galaxy stellar mass function (GSMF) and 
halo mass function measured in simulations: 

n g (>M starii ) = « A (>M vir)i ) (1) 

where n g and «/, are the number density of galaxies and ha- 
los (note again that "halos" here and throughout refers to both 
distinct halos and their subhalos), respectively. Galaxy stel- 
lar masses can hence be assigned to halos at any epoch once 
the GSMF at that epoch is known. In the simplest version of 
this scheme there is assumed to be no scatter in the relation 
between halo mass and stellar mass (see i j4.5| for a discussion 
of scatter in the context of our models). 

Subhalos lose mass due to tidal stripping as they orbit 
within their parent halo. Since stripping primarily affects the 
outer regions of the subhalo, we expect the galaxy, which re- 
sides within the inner few kpc of the subhalo, to be relatively 
unaffected by this process. Hence, for subhalos, when relat- 
ing halo mass to luminosity or stellar mass we use its virial 
mass at the epoch when it is first accreted onto the parent 
halo, MjJ? c , rather than its mass at the epoch of observation. 
This ch oice is well motivated b oth by hydrodynamical simu- 
lations (Nagai & Kravtsovl2005l) and detailed modeling of the 
small to intermediate scale (0.1 < r < 10 h~ l Mpc) clus tering 
of galaxies over a range of redshifts dConrov et alj|2006h . 

One may ask to what extent it is justifiable to identify satel- 
lite galaxies with subhalos in dissipationless simulations. It 
has been shown that the subhalo population in dissipation- 
less simulations is indeed quite similar to the galaxy pop- 
ulation in hydrodyn amical simulations (Zheng et alj l2005t 
Weinberg et al.l l2006) and semi-analytic models dzheng et all 
20051) . In particular, satellite populations in these hydrody- 
namical simulations with cooling and galaxy formation have 
an almost identical halo occupation distribution to the subha- 
los in dissipationless simulations. These conclusions are cor- 



roborated by the general success of the subhalo-based models 
of galaxy clustering. 

This model, for example, accurately captures observed re- 
lations between cluster luminosity and the nu mber of galaxies 
within a clu ster as a function of cluster mass (Vale & Ostriker 
12004 120061) . the luminosity dependence of the galaxy-matter 
cross-correlation function (after a reas onable amount of scat- 
ter is introduced into this relat ion; see [j asitsiomi et al] |2004 
for details), close pair counts dBerrier et al.ll2006l) . the lumi- 
nosity, scale, and redshift-dependen ce of the galaxy au tocor- 
relation function from z ~ 5 to z ~ dConrov et al.ll2006|). a nd 
mass-to-light ratios in local clusters (Tasitsi omi et alj |2007). 

Emboldened by the success of this simple model, in the 
present work we extend it by populating halos with galaxies 
at z ~ 1 in the above way and then using halo merger trees 
derived from A^-body simulations to follow the evolution of 
these galaxies to z ~ 0. Specifically, we use the observed 
galaxy stellar mass function (GSMF) of iFontana et al.l (2006) 
at z ~ 1 to assign stellar masses to halos, and then follow 
the dissipationless evolution of these galaxies via the merg- 
ing history of their dark matter halos to z ~ 0. This exercise 
is appropriate for the evolution of the most massive galaxies, 
as these galaxies formed the bulk of their stars at z > 2 (e.g. 
[ Bower et al.l 1 1 992h iTrager et al.ll2000l Ivan Dokkum & Franxl 
I200U iThomas et alj 120051: Uimenez et all 120061) and hence 
largely evolve dissipationlessly at z < 1. Note that neglect- 
ing star formation at z < 1 means that the evolution of the 
total stellar mass predicted by our models is a lower limit on 
the actual evolution. 

The follo wing z = 1 GSMF Schechter parameters are 
adopted from Fonta naet alj (|2P06): 

a = -1.26±0.1 

M* = 10 1101±01 M Q (2) 
0* = (7.6 ± 2.4) x 10- 4 Mpc" 3 . 

The la err ors quoted above ar e not the statistical errors 
reported in Fontana et al] ( 120061) but are instead meant to 
roughly encompass the various published estimates of the 
z = 1 GSMF Schechter parameters. The statistical errors are 
a factor of 2 - 3 smaller than these approximate systematic 
errors. Here and throughout the Chabrier IMF is used when 
quoting stellar masses. 

It is important to note that while different authors derive 
somewhat different Schechter parameters for the z ~ 1 GSMF, 
all authors agree that the massive end of the GSMF (> M„) 
evolves very little, if at all, since z ~ 1. In addition, while 
several measurements of th e GSMF at z ~ 1 h ave relied on 
photometric redshifts (e.g., iBorch etaT] 120061) . the general 
conclusions from these studies have been support ed by mea- 
surements which utilize sp ectroscopic redshifts dBundv et alj 
20061: IFontana et alj|2006h . 

Figure Q] presents a comparison betwee n various observe d 
GSMFs at z ~ 1 and the z - GSMF from lCole et aL1(l200lh. 
Note that the z ~ 1 GSMF used in our models (IFontana et alj 
120061) is below all the other z - 1 GSMFs. This implies that 
our results concerning the z = brightest cluster galaxy (BCG) 
luminosities derived with this GSMF will be lower bounds 
relative to the other GSMFs. 

At z = 1 , a space density of 1 0~ 4 Mpc" 3 corresponds to halos 
of virial mass M sn ~4x 10 13 M Q and stellar mass M star ~3x 
10 u Mq. At this space density, ~ 17% of halos are subhalos 
at z ~ 0. Finally, note that although we formally assign stellar 
mass to all the halos found in the simulations and track their 
evolution as described in the following section, our results are 
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FIG. 1. — Observed galaxy stellar mass functions (GSMFs). The non- 
solid lines are best fit Schechter functions to the following GSMFs at z ~ 
1: IBorchetalJ KOOft dot-dashed line), IBundv et all <2006L do t ted lin e), 
IFontana et alj (2006L dashed line). The GSMF from lProrv et al.l 120041) is 
plotted directly, both for their fiducial GSMF (diamonds) and an estimate 
which includes lost-light corrections (solid circle s; error bars a re Poisson un- 
certainties only). The GSMF at z ~ from Cole et al. 12001, solid line) is 
also included, along with arrows indicating a change of 0.2dex in M st ar at 

= 10- 4 . 

quite insensitive to low mass halos and are instead governed 
by the evolution and fate of much more massive halos. Since 
it is these halos that are most easily resolved and tracked from 
timestep to timestep, we expect our results to be insensitive to 
our simulation resolution. 

3.2. Dynamical evolution models 

Our model contains (at least) three adjustable components. 
The first is our assumed cosmology and will be fixed through- 
out, though we comment qualitatively on the effects of chang- 
ing certain cosmological parameters in i!4.5l The second com- 
ponent is our input method for assigning stellar masses to ha- 
los at z ~ 1 ■ In Sj43]we demonstrate that our results are robust 
to reasonable changes of this second component. In partic- 
ular, we introduce scatter in the assignment between stellar 
mass and halo mass and marginalize over the uncertainties in 
the z ~ 1 GSMF and find no qualitative change to our conclu- 
sions, provided that the true GSMF is not considerably differ- 
ent from recent estimates. 

The third adjustable component is the most uncertain and 
concerns how we treat subhalos that have dropped out of the 
halo catalog. This can occur either because the subhalo is 
physically disrupted or because it is stripped below the reso- 
lution limit of the simulation; it is often quite difficult to dis- 
tinguish these two cases within the simulation. We do nothing 
to the galaxies assigned to subhalos while the subhalos re- 
main identifiable in our simulation — i.e., the satellite galaxy 
within the subhalo experiences no tidal stripping. Once a sub- 
halo is destroyed, we are free to redistribute the stars from 
this subhalo in one or more of the ways outlined in the In- 
troduction. To summarize that discussion: the stars can be 
deposited onto the central galaxy (this assumes that the satel- 
lite within the destroyed subhalo has merged with the central 
galaxy), can remain as a satellite galaxy without an identifi- 
able subhalo, and/or the stars can be added to the intra-cluster 
light (ICL). In the latter case the stars are added to the outer 
regions of the central galaxy, beyond the optical radius (the ra- 
dius within which the central galaxy luminosity is measured). 

In order to explore these possibilities, four models are con- 
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Summary of Models 


Model 


Fate of satellite galaxy 




in a disrupted subhalo 


Sat2Cen 


Stars deposited onto the central galaxy 


KeepSat 


Stars remain bound as a satellite galaxy 


Sat2ICL 


Stars deposited into the ICL 


Sat2Cen+ICL 


Stars divided equally between the ICL 




and the central galaxy 



structed which differ only in the fate of the stars within dis- 
rupted subhalos. Model Sat2Cen assumes that all of the stars 
are deposited onto the central galaxy. Conversely, model 
KeepSat assumes that the stars remain bound as a satellite 
galaxy. Model Sat2ICL assumes that all stars are deposited 
into the ICL. Finally, model Sat2Cen+ICL assumes that the 
stars are distributed equally between the ICL and the central 
galaxy. We will assume below, when comparing to observa- 
tions of the ICL, that the ICL in these models is generated 
predominantly from the remnants of mergers with the cen- 
tral galaxy at z < 1. These assumptio ns are motivated by 
hydrodynamical simulations of clusters ( Willman et al| 2004 ; 
[ Murante et aTll2004t iRudick et alJl2006tlMurante et alj|2007 : 
Somm ~Larsen etal.1 12005b which showed that the majority 
of the ICL is built up at z < 1 from major mergers with the 
central galaxy, rather than tidal stripping as the satellite orbits 
in the cluster. For reference, these models are summarized in 
Table Q] 

These models have implicitly assumed that all stellar mass 
at z = 1 persists to z = 0. However, stellar mass-loss due to 
winds and supernovae can result in a significant decrease in 
the aggregate stellar mass of a population over time. In or- 
der to understand these effects, consider a secular mass-loss 
rate of M/M = 0.05(f/Gyr) _I for a stellar population formed 
in an instantaneous b urst and older than a few hundred Myr 
dJungwiert et alj|2001l) . If massive red galaxies grew most of 
their stellar mass in the form of a single burst at z = 2, then 
the fraction of stellar mass lost between z = 1 and z = is only 
7%; however, if the stars in these galaxies all formed at z = 1 
then the fraction is 36%. Since observations place the epoch 
of star-formation in these massive galaxies at z > 2, mass-loss 
effects are likely unimportant. 

3.3. Generating luminosities 

Comparison to observations in $4.2\ and ^4.31 will require 
conversion from stellar masses to A'-band and /-band lumi- 
nosities, respectively. To make this conversion we use the 
relation between mass-to-light ratios and colors provided by 
iBell etalJd2003h . 

M stal /L K = 0.72 

M stm -/L, = 1.90 (i> 

where we have assumed a color of (u - r) = 2.5 and (B-R) = 
1.5 when deriving M s1m /Lk and M stal /Lj, respectively, and a 
Chabrier IMF. These colors are appropriate for the bright end 
of the red sequence dPahrdl999HBaldrv et al.l2004t[Bell et alj 
120031) . 

All the results to be discussed below are independent of 
our assumed IMF because we are only interested in relative 
evolution from z = 1 to z = 0. Thus, so long as the observa- 
tions at these epochs use the same IMF, the results are insensi- 
tive to the particular IMF used (whether for example Chabrier, 
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Kroupa, or Salpeter IMFs are used). Results which concern 
luminosities are also IMF independent so long as we use a 
mass-to-light ratio with the same IMF as used in our GSMFs 
(as we have done above). 

The luminosities of galaxies at z ~ in these models are 
strictly lower bounds for two reasons: 1) these models ne- 
glect star-formation since z = 1, and 2) some galaxies might 
not be as red as the colors assumed in the previous paragraph, 
and hence they will be brighter at a fixed stellar mass. For 
the massive BCGs studied herein however, these effects are 
unimportant. It will become apparent in the next section that 
including these possibilities would only strengthen our gen- 
eral conclusions since both residual star-formation and bluer 
colors would increase the BCG luminosities and the evolution 
of the massive end of the GSMF 

We also assume that the ICL has the same color as the tip 
of the red sequence when converting ICL stellar mass to lumi- 
nosities. Such an assumption appears borne out by observa- 
tions of an at most weak color gradient out to several hundred 
kpc fro m the BCG, with some authors finding a slight red- 
dening dGonzalez et alj |2000; Krick et al. 2006) a nd others a 
slightly bluer color, or no gradient at all "dZibetti et alj |2005). 

4. RESULTS 

We now compare the models constructed in |J3] to the ob- 
served evolution in the galaxy stellar mass function ( 34. U . to 
the relation between BCG luminosity and cluster virial mass 
at z ~ ( £|4.2| >, and to properties of the intracluster light (ICL; 
H4.3\ . In i!4.4l we discuss these models in the context of star 
formation since z = 1 ; at the end of this section several caveats 
and assumptions made herein are explored. 

4.1. Evolution in the GSMF 

In this section we present the evolution in the galaxy stellar 
mass function (GSMF) for the models described in £j3j and 
compare to observations. An effective way to quantify the 
evolution in the massive end of the GSMF is by quantifying 
the change of stellar mass corresponding to galaxies with a 
fixed given value of spatial number density. Hence, we define 
the quantity M_4, such that <f>(M-4) = 10~ 4 , and its evolution 
as AM_ 4 = APf-Mf^ 1 (see the top panel of Figure [2] for a 
pictorial representation). This quantity is more stable than its 
inverse (the evolution in the number density of a given stellar 
mass) due to the exponential cut-off of the GSMF at high stel- 
lar masses. Note that M-4 is dominated by the space density 
of massive, rare objects. 

Uncertainties in the observed z ~ 1 GSMF are incorporated 
into the models by generating 200 realizations of the GSMF 
with each Schechter parameter drawn from a Gaussian distri- 
bution with mean and dispersion equal to the best fit and la- 
errors on the observed GSMF (see Equation[2]l. We have tried 
other GSMFs at z ~ 1 dBundv et alJl2(M iBorch et al.ll200^ 
and find qualitatively similar results. 

The bottom panel of Figure [2] shows the distribution of 
AM_4 generated from the 200 Monte-Carlo realizations for 
models Sat2Cen and Sat2Cen+ICL (models KeepSat and 
Sat2ICL produce no/little change in AM-4, see below). The 
observational limits are denoted by the shaded band and have 
been estimated from observations of the z = 1 GSMF (see Fig- 
ure [TJ. Note that for each model the predicted evolution is a 
lower bound, since star-formation between z ~ 1 and z ~ 0, 
which is neglected in these models, will increase AM-4. How- 
ever, since the massive galaxies which dominate this quantity 
have formed the bulk of their stars at z > 1, we expect the 
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FIG. 2. — Evolution in the GSMF from z ~ 1 to z ~ 0. Top panel: Pic- 
torial representation of how we quantify evolution in the massive end of the 
GSMF with the parameter AM_4, for one realization of model Sat2Cen (error 
bars denote Poisson uncertainty only). See the text for details. Bottom panel: 
Distributions of AM_4 for models Sat2Cen and Sat2Cen+ICL . These distri- 
butions are obtained by running each model using a slightly different GSMF 
to populate halos at z ~ 1 ■ Model KeepSat predicts no change in the entire 
GSMF and hence AM4 = for that model, while model SatHCL results in 
AM_4 ~ 0. Running models for a series of GSMFs demonstrates the effect 
of uncertainties in the z ~ 1 GSMF on our results. Observational limits are 
denoted by the shaded region. 



contribution to AM-4 from star-formation to be unimportant. 
Note that the evolution in the GSMF presented in Figure [2] 
makes no reference to the observed z ~ GSMF. We simply 
compare the observed GSMF at z = 1 with the GSMF evolved 
to z = with our models. For comparison, AM-4 = 0.4 for 
dark matter halos in our adopted cosmology. 

For each model, the differences in AM-4 for different 
Z ~ 1 GSMF Schechter parameters arise because different 
Schechter parameters result in a different relation between 
stellar mass and dark matter halo mass, via Equation [T] The 
growth of a dark matter halo and, through our models, the 
growth of the central galaxy, is driven by accreted halos span- 
ning a range in mass. Therefore, the predicted growth of the 
central galaxy mass will change, if the mapping between halo 
mass and stellar mass changes. We now explain the behavior 
of each model in turn. 

Model Sat2Cen displays the largest increase in M-4 be- 
cause massive galaxies, which by definition dominate this 
quantity, are growing rapidly. In this model, galaxies with 
Mstar > 1O U M0 at z = have on average more than doubled in 
mass since z = 1. Rapid growth of massive galaxies occurs in 
this case because the satellite galaxies within disrupted sub- 
halos add all of their stars onto the central massive galaxy. 
Model KeepSat predicts AM-4 = 0, since in this model galax- 
ies do not evolve, i.e., galaxies neither merge nor form stars 
since z = 1. In Model Sat2ICL the massive end of the GSMF 
does not increase (but in some realizations actually decreases 
slightly with time as some very massive satellite galaxies dis- 
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FIG. 3. — Lumino sity of the BCG as a function of the cluster virial mass, 
comparing data from Lin & Mohr (2004) (solid circles) to models (open dia- 
monds). 



rupt) because the stars within satellite galaxies are transferee! 
to the ICL, which, for the purposes of the GSMF amounts to 
deleting the galaxy from the sample. Model Sat2Cen+ICL is, 
by construction, intermediate between models Sat2Cen and 
Sat2ICL, since half of the stars from disrupted subhalos are 
deposited onto the central galaxy and the rest into the ICL. In 
this model the most massive galaxies have increased in mass 
by the more modest factor of ~ 50% since z = 1 . 

Based on the evolution in the GSMF, model Sat2Cen 
is strongly disfavored. Models KeepSat, Sat2ICL, and 
Sat2Cen+ICL fair far better. In fact, based on current ob- 
servations, all of these models appear more or less equally 
viable (their relative viability depends on what one assumes 
about the observationally allowed range in AM-4). We now 
turn to comparisons of the models with observations of BCGs 
and the ICL, in the hope of more strongly distinguishing be- 
tween them. 

4.2. BCG Luminosities at z ~ 

We now confront our models with obser vations of clus- 
ter properties at z ~ 0. iLin & Mohrl (|2004) have computed 
virial masses from X-ray observations and, using the 2MASS 
database, estimated BCG luminosities 5 for 93 clusters at z < 
0.1. Figure Opiots the /T-band BCG lumino s ity ver sus cluster 
virial mass for the data from ILin & Mo hr (2004, solid cir- 
cles) and for the models 6 . Several trends are apparent. Mod- 
els KeepSat and Sat2ICL predict identical BCG luminosities 
(since in both models the central galaxy, which is identified as 
the BCG, does not accrete any satellite stars since z = 1) and 
are in good agreement with the observations. Model Sat2Cen 
is, again, in strong disagreement with the observations. Fi- 
nally, the predictions of model Sat2Cen+ICL are in between 

5 Note that ~ 70% of the BCGs in this sample are centered in the cluster 
to within 5% of the virial radius; in what follows we assume that all BCGs in 
this sample are the central galaxy. 

6 For this comparison, we have converted our halo mass defi nition of 334 
times the mean density of the Universe to the definition used in Lin & Mohr 
12004, 200 times the critical density), by using an NFW density profile with 
a mass-dependent concentration. 
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FIG. 4. — Top panel: Absolute /-band magnitudes of the combined BCG 
and ICL components as a fu nction of total cluster mass. Observations at z ~ 
(Gonzalez et al. 2005, 2007, solid circles) are compared to models KeepSat 
(squares), Sat2Cen, Sat2ICL, and Sat2Cen+ICL (crosses). Bottom panel: 
fraction of BCG+ICL luminosity that is contained in the ICL component as 
a function of cluster virial mass, comparing observations to models Sat2ICL 
(diamonds) and Sat2Cen+ICL (triangles). Models Sat2Cen and KeepSat pro- 
duce no ICL component. 



those of models KeepSat and Sat2ICL and model Sat2Cen, by 
construction, and are mildly disfavored by the observations. 

The failure of mod el Sa t2Cen is of course no sur prise in 
light of the results in E l4-.ll The failure is, as in &I4.2I simply 
a manifestation of the fact that the massive end of the halo 
mass function in a ACDM cosmology evolves much more 
strongly from z = 1 to z = than the observed evolution in 
the GSMF. This is corroborated by observational constraints 
on halo masses at various epochs which indicate that while 
the stellar and dark mat ter components grow in lock-step 
for lo wer-mass systems dHeymans et alJl2006t IConroy et alJ 
120071) . the stellar mass growth of central galaxies in high mass 
halos appe ars to be outpaced by the growth in their halo mass 
at z < 1 dConrov et aHl2007l) . 

Comparison to BCG luminosities does however provide 
stronger constraints on model Sat2Cen+ICL compared to the 
constraints from evolution in the GSMF. Specifically, the dis- 
agreement between observations apparent in Figure [3] sug- 
gests that less than half of the stars from disrupted subhalos 
can end up in the central BCG. And indeed, from the agree- 
ment between observations and models KeepSat and Sat2ICL, 
we conclude that no growth is favored. 

4.3. The ICL Component at z ~ 

We now confront our models with observations of the ICL 
component. The notion of intracluster light arose from the 
observation that the extend ed profiles of BCGs were in excess 
of de Vaucouleurs profiles (Matthews et al. 1964; Schombert 
1988). There is currently no strong consensus on whether the 
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ICL is simply the outer component of the BCG or whether it 
is dynamically distinct. Observationally, the ICL is often de- 
fined as the total light beyond a particular surface brightness 
level, although recently there have been attempts to model the 
entire surface brightness profile with multiple components, 
thus separating the BCG and ICL in a less arbitrary way. 
For o ur purposes we use the data from lGonzalez et ail d2005l 
2007) who have measured the surface brightness profiles for 
24 BCGs at z < 0.12 in the /-band and have also measured 
virial masses for the clusters. 

When considering the ICL, the most straightforward ob- 
servable to confront with models is the combined BCG and 
ICL light. For this quantity one does not need to rely 
on the potentially arbitrary distinction between BCG and 
ICL. The top panel of Figure |4] plots the absolute /-band 
magnitude for the combined BCG and ICL components, 
Mbc g+icl, as a funct i on of cluster virial mass 7 , for the data 
from iGonzalez et al.l d2005l solid circles), and for our mod- 
els. Models Sat2Cen, Sat2ICL, and Sat2Cen+ICL predict the 
same Mbcg+icl, since these models differ only in the way in 
which the stars are distributed between the BCG and ICL. As 
can be seen from the figure, these three models all adequately 
reproduce the observations over a range of cluster masses. 

Model KeepSat however predicts a substantially different 
A^bcg+icl, since in this model no stars are added to the ICL 
nor BCG since z = 1. In particular, model KeepSat pre- 
dicts Mbcg+icl > lmag lower than observations, which cor- 
responds to a factor > 2.5 lower luminosity, and because of 
our adopted constant mass-to-light ratio, this corresponds to 
the same factor lower in stellar mass. This discrepancy is far 
too great to be accounted for by the small effects neglected 
in these models, such as star-formation, tidal stripping, and 
ICL generation at z > 1 (see below). Hence, observations of 
A^bcg+icl strongly suggest that model KeepSat is unrealistic. 

A more complex, but potentially more discriminating ob- 
servable, is the fraction of BCG and ICL light that is in the 
ICL. In this case comparison between models and data must 
be treated carefully because the separation between ICL and 
BCG is not handled in the same way for different datasets. 
Our operational definiti on of ICL i s sim ply the light not 
counted as the BCG by lLin &Mohri d2004l) . Their definition 
of BCG luminosity corresponds to the light within a surface 
brightness of fix ~ 21 mag arcsec" 2 . Assuming I-K = 2, 
which is appropriate for bright red galaxies, implies a sepa- 
ration between ICL and B CG at /j,/ = 23 mag ar csec" 2 . The 
observational results from Gonzal ez et alj d2005i) . which are 
in the /-band, have been recast in this way to afford the most 
robust comparison to our model (A. Gonzalez, private com- 
munication). 

The bottom panel of Figure H] plots the ICL fraction 
for models SatllCL and Sat2Cen+ICL (recall that models 
Sat2Cen and KeepSat do n ot have an ICL comp onent) and 
compares to the results from Gonz alez et al.l d2005l solid cir- 
cles). It is clear that model Sat2ICL predicts much more ICL 
light than the other models and an ICL light fraction that is in 
excellent agreement with the observations. 

There are two additional routes by which the ICL can be 
built up that have been neglected thus far: build up of the ICL 
at z > 1 and the tidal stripping of satellites as they orbit within 

7 Their definition of virial mass is the mass enclosing a region with mean 
density equal to 500 times the critical density; we have converted bo th th eir 
masses and our to a definition of 200 times the critical density; see ^4.2l for 
details. 



the cluster potential. Specifically, hydrodynamical simula- 
tions have found that ~ 85% of the stars in the ICL at z = 
were deposited at z < 1 ( Wil lman et al.ll2004t iMurante et al.1 
120071) . and less than 30% of the ICL was built u p by tidal strip- 
ping of satellite galaxies (Murante et al. 2007). The majority 
of the ICL is thus built during violent merging events with 
the central galaxy and/or the complete disruption of satel- 
lites — i.e. the two processes that are captured in our treat- 
ment of the ICL. Recent observations of the color of the ICL 
suppo rt the picture that stars comprising the ICL formed at 
Z > 1 dKrick et aUl2006l) . While both processes will increase 
the ICL component, neither will produce enough additional 
ICL to account for the discrepancy between models Sat2Cen, 
KeepSat and Sat2Cen+ICL and observations depicted in Fig- 
ure [4] — if these simulations are accurately capturing the 
build-up of the ICL. 

4.4. Implications for Star Formation Since z = 1 

Until now we have focused on observations that can be de- 
scribed with purely dissipationless modeling. Now that we 
have identified a dissipationless model that adequately repro- 
duces various observations (model Sat2ICL), we can ask what 
more must be added to such a model in order to reproduce the 
observed global galaxy population. Clearly, the most impor- 
tant process neglected thus far is star formation, which be- 
comes increasingly important in lower mass halos. We now 
turn to a discussion of the importance of star formation as a 
function of halo mass since z = 1 . 

We first construct a simple model that places the "true" z = 
stellar mass in dark matter halos. This is accomplished us- 
ing the methodology outlined in |j3j now matching the z = 
GSMF to the z = halo mass function. Such a model will 
have the correct z = GSMF by construction and should have 
approximately the correct relation between stellar mass and 
halo mass since this method has been shown to reproduce nu- 
merous observations remarkably well (see fj3]for details). 

Figure|5]compares these true stellar masses to stellar masses 
from models Sat2Cen, Sat2ICL and Sat2Cen+ICL as a func- 
tion of z = halo mass (top panel) and stellar mass (bottom 
panel). As before, we generate 200 realizations of these mod- 
els by sampling the z = 1 GSMF uncertainties (recall that in 
these models the z = stellar masses are products of the z = 1 
GSMF combined with the dark matter halo merger trees to 
z = 0). The resulting mean and la dispersions are included in 
this figure. At large masses, the stellar masses from models 
Sat2ICL and Sat2Cen+ICL match the "true" stellar masses 
while model Sat2Cen overpredicts the true stellar masses, 
although all are consistent with the true masses at roughly 
the 2er level. This is not surprising, both in light of the re- 
sults from previous sections and more generally because mas- 
sive galaxies (which reside in massive halos) are observed 
to have finished forming stars by z ~ 2. In fact, if there is 
truly zero star-formation in these massive galaxies since z = 1, 
then Figure [5] suggests that the z = "true" stellar masses are 
best reproduced by a model in between models Sat2ICL and 
Sat2Cen+ICL, i.e. subhalos transfer perhaps ~ 80% of their 
stars to the ICL, and ~ 20% to the central galaxy. A distinc- 
tion this refined should, of course, be treated with caution. 

The behavior at lower masses is more interesting. In this 
regime the stellar masses from models Sat2Cen, Sat2ICL and 
Sat2Cen+ICL are substantially less than the true masses, indi- 
cating that the z < 1 star formation is increasingly important 
in halos of lower mass. In fact, according to these models 
roughly 40% of the stars in z = halos of mass M v ; r ~ 10 13 
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FIG. 5. — Ratio of the z = stellar masses predicted by models Sat2Cen, 
SatHCL and Sat2Cen+ICL to the "true" stellar masses, as a function of z = 
host halo mass (top panel) and stellar mass (bottom panel). These "true" 
stellar masses are obtained by matching the observed z = GSMF to the 
Z = halo mass function (see the text for details). Dotted lines denote the 
lcr dispersion around model Sat2Cen+ICL estimated from 200 Monte Carlo 
realizations that incorporate the z = 1 GSMF uncertainties (the dispersion 
around the other models is similar). 

h _1 M Q formed at z < 1 while in halos of mass M v j r ~ 10 12 
h _1 M Q ~ 80% of the stars were formed over the same inter- 
val. This is general l y consistent with o bserved trends (e.g., 
iHeavens et all 12004 iNoeske et al.ll2007l) . The implied star- 
formation rates from this comparison are sensitive to the less 
massive end of the GSMF at z = 1, which is less-well con- 
strained observationally, and should thus be treated with a 
caution. Note that the models converge at lower masses, and 
thus these conclusions are insensitive to the way in which dis- 
rupted subhalos are handled. 

4.5. Caveats & Assumptions 

4.5.1. Definition of the BCG 

It has recently come to light that standard photometry in 
large galaxy surveys systematically under estimates BCG lu- 
minosities for several different reasons (lLauer et alj |2007; 
iBernardi et af]|2007b . This issue is complicated by the some- 
what arbitrary distinction between BCG and ICL, as one 
needs a well defined notion of a BCG in order to claim that 
standard techniques are "missing" BCG light. The effect can 
be as large as 1 magnitude, though in such cases it appears 
that the ICL is included as part of the BCG. Our results are 
robust against these effects because our notion of a BCG is 
precise ly that measured b y the data to which we compare 
(i.e. bv lLin & Mohril2004l) . while the ICL is simply light out- 
side the optical radius (i.e. outside the region counted as the 
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FIG. 6. — Stellar mass as a function of virial mass at z = 1. The solid 
line indicates our standard one-to-one correspondence between halo and stel- 
lar mass, while the symbols represent a prescription for scatter between the 
two quantities (only 25% of the total number of objects are plotted for clar- 
ity). The dashed and dot-dashed lines represent the stellar mass that a galaxy 
would have if it had a star-formation efficiency of r] = 1.0 and 0.1, respec- 
tively (assuming that the baryon-to-dark matter ratio in each halo is the cos- 
mic mean value of 0.17). See the text for details. 

BCG). Unfortunately, this means that our results which rely 
on the separation between the BCG and ICL are not directly 
exportable to other observations of the BCG and ICL if such 
observations separate these two observables in different ways. 
As mentioned previously, a more robust approach to this type 
of modeling would be to present actual surface brightness pro- 
files which could then be compared to any we ll-defined obser- 
vational sample (see e.g jRudick et alj[2006l) . 

4.5.2. Scatter in the galaxy-halo connection 

We now explore the impact of scatter in the M^-M^ re- 
lation on our results (recall that the relation between stellar 
mass and halo mass utilized in the previous section was gen- 
erated by assuming a one-to-one correspondence). In order 
to explore the maximal effect that scatter can have on our 
results, we generate a rather extreme prescription of scatter. 
The scatter is included by multiplying each halo mass by a 
random number drawn from a Gaussian with a = 0.6dex (see 
Tasitsiomi et al. 2004], for a different prescription of scatter). 
Stellar masses are then matched to this random number in the 
standard way 8 . Figure|6]compares this model (solid symbols) 
to our standard one-to-one correspondence (solid line). This 
figure includes both distinct halos and subhalos; for the latter 
we use the virial mass at the epoch of accretion, as before. 
In addition, we include lines that indicate the amount of stel- 
lar mass a galaxy would have if its halo contained the cosmic 
mean baryon-to-dark matter ratio, ft, = 0.17, and it converted 
a fraction rj of those baryons into stars (rj is often called the 
star-formation efficiency). 

Comparison between the scatter prescription and the lines 
of constant star-formation efficiency indicate why this pre- 
scription is extreme — there are galaxies which have rj ~ 1 
and indeed some rare cases where the baryon fraction in the 
halo exceeds ft,. At z ~ the star-formation effic iency is 
almost certainly rj < 0.25 (Mandelbau m et aDl2006l and de- 
creases for increasing halo mass), so we have effectively scat- 

8 If £ denotes the halo mass multiplied by a random number drawn from 
a Gaussian, then the stellar masses are assigned via n g (> M star \) = «;,(> 
c.f. Equation[T] 
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tered galaxies, at a fixed stellar mass, to halo masses that vi- 
olate, or almost violate constraints on both 77 and //,. Note 
that this type of scatter (scattering down in halo mass at a 
fixed stellar mass) will likely cause the most change in the ob- 
servables discussed in the previous sections because massive 
galaxies in less massive halos will have less violent accretion 
histories since z = 1 compared to more massive halos. 

However, even with this large amount of scatter, our re- 
sults remain qualitatively unchanged, though our conclusions 
are not as strong. For example, the model predictions for 
the LftCG-M vn relation are lower by ~ O.ldex, with a more 
pronounced tail toward lower Lbcg at lower masses. Evo- 
lution in the GSMF is also less by about the same amount. 
The BCG+ICL luminosity-mass relation still rules out model 
KeepSat, and the ICL fraction still favors model Sat2ICL. 

4.5.3. Observational Uncertainties 

In £14.1 1 we explored the sensitivity of the model predictions 
to the adopted z ~ 1 GSMF; we now do the same for the BCG- 
cluster mass relation presented in S14.2I We have checked 
by eye the BCG-cluster mass relation produced from model 
Sat2Cen for the 200 realizations of the z ~ 1 GSMF a nd com- 
pared them to the observations of iLin & Mohrl (120041) . Model 
Sat2Cen matches the observations approximately ~ 13% of 
the time (note that this fraction is the same fraction of real- 
izations of model Sat2Cen that match the observational con- 
straints on evolution in the GSMF, see 34.11 1, In these cases 
either M» or 0* (or both) are > lcr below the mean GSMF 
Schechter parameters. Recall that these uncertainties are 
rough combined systematic and statistical uncertainties and 
are th us likely upper bou nds. Further, the GSMF we adopt 
from iFontana et al. (2006) is a lower bound with respect to 
other observations at z ~ 1 (see Figure Q}. Nevertheless, any 
future revision of the GSMF downward by > 1 a (systematic + 
statistical) would weaken the constraints on the models signif- 
icantly. Clearly, our results and conclusions could be sharp- 
ened with more accurate measurements of the GSMF in the 
future. 

4.5.4. Effects of cosmology 

In our simulations of the ACDM cosmology, the adopted 
normalization of the power spectrum, a% = 0.9, was some- 
what higher than rece nt constraints from the 3-year WMAP 
data: cr 8 « 0.75-0.8 (ISpergel et al.ll2006l) . A lower value of 
erg would imply that the same observed galaxy number density 
(and stellar mass) corresponds to smaller halo circular veloc- 
ity and virial mass. There would also be more evolution in the 
abundance of massive halos between z = 1 and z = 0, as halo 
formation times would be shifted to more recent time. This 
would imply that the amount of evolution would be larger 
and our conclusions would be stronger in a lower erg universe. 
However, later formation times also imply that there would 
be less time for accreting halos to merge and contribute to the 
growth of the BCG stellar mass. The relative importance of 
these competing effects will have to be quantified in a future 
analysis similar to the one presented here with simulations of 
lower erg. 

5. DISCUSSION 

5.1. Implications 

We have explored four models for the dissipationless evo- 
lution of galaxies since z = 1 . These models were constructed 
to match the galaxy stellar mass function (GSMF) at z = 1 and 



differ only in the fate of satellite galaxies when the subhalo, 
within which the satellite is embedded, is disrupted. Upon 
confronting these models with various observations we have 
found that only a model in which a significant fraction of stars 
(> 80%) from disrupted subhalos are transfered to the ICL 
(referred to as model Sat2ICL above) is consistent with data. 
The failure of the other models provides significant insight 
into the dissipationless evolution of galaxies. 

A model in which all the stars from disrupted subhalos are 
transfered to the central galaxy (model Sat2Ceri) is strongly 
ruled out both by observations of the evolution in the GSMF 
and observations of the z = Lbcg _ A^vir relation. Such a 
model would only be viable if the observed GSMF at z ~ 1 
were significantly revised downward from the current low- 
est reported measurements. The failure of such a model im- 
plies that, if stars from disrupted subhalos are transfered to the 
central galaxy, then they cannot be put in the central regions 
where BCG luminosities are measured. Such a conclusion is 
corroborated by simulations of dissipationless ("dry") galaxy- 
galaxy mergers, which find that the resulting galaxy gener- 
ally becomes mor e extended rather than subst antially brighter 
at the center (e.g. iBovlan-Kolchin & Mall2007l although this 
conclusion depends on the orbital parameters of the accreted 
satellite). 

We have also explored an extreme model where satellites 
never disrupt even when their subhalos do (model KeepSat). 
This model fails dramatically when compared to observations 
of the combined luminosity of the BCG and ICL, under the 
assumption that the ICL is built-up predominantly at z < 1 . 
Although this assumption appears justified in light of recent 
hydrodynamical simulations, one should note that our con- 
clusion as regards model KeepSat relies on this assumption. 
In our models the massive subhalos correspond to massive 
satellites and it is the fate of these massive satellites that most 
strongly affects the comparison to observations of the com- 
bined luminosities of the BCG and ICL. The failure of this 
model thus strongly suggests that the disruption of subhalos 
in our high-resolution iV-body simulations corresponds to the 
disruption of satellite galaxies, at least for the most massive 
subhalos. This need not have been the case; recent semi- 
analytic models (SAMs) have decoupled the dynamical evo- 
lution of subhalos from satellites wh en the subhalo disrupts 
dCroton et al.ll2006l;IWang et al.ll2006l) . and hence these mod- 
els produce a significant population of satellites with no iden- 
tifiable subhalo (the so-called "orphan" population). 

The failure of model KeepSat strongly suggests that in fact 
satellites disrupt when their subhalo disrupts, at least for mas- 
sive satellites, and hence any model (including the SAMs 
mentioned above) which fails to tie the fates of massive satel- 
lites to their subhalos will likely fail to reproduce the ob- 
served combined luminosities of the BCG and ICL. 9 More- 
over, the failure of model KeepSat provides additional justifi- 
cation for previous modeling where a one-to-one relation be- 
tween galaxies and hal os extracted from jV-body simulations 
has been assumed (e.g.lKravtsov et al.l l2004; Tasit siomi et alj 
l200ltlConrovetalJl2006l) . " 

The data instead favor models where most, if not all, of 
the stars from disrupted satellites are deposited into the ICL. 
Model Sat2ICL puts all stars from disrupted satellites into the 

9 However, the simulations used in the current generation of SAMs have 
coarser resolution than the simulations used herein. It is thus difficult to make 
a fair comparison between the different approaches to handling the satellite 
— subhalo correspondence. 
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ICL while model Sat2Cen+ICL puts only half into the ICL 
and the rest onto the BCG. Although comparisons with vari- 
ous observations favor model Sat2ICL, the uncertainties and 
assumptions discussed in previous sections suggest that re- 
ality may lie somewhere in between these two models. The 
strongest discriminant between these scenarios is the fraction 
of BCG+ICL light contained in the ICL; as elsewhere, model 
Sat2ICL most faithfully reproduces these observations. 

There is a growing consensus that ma ssive red galaxies 
were more or less in place by z ~ 1 (e.g. IWake et all 120061: 
ICimatti et aLll2006t iBundv et ai1l2006l) . At first glance it ap- 
pears difficult to reconcile this fact with ACDM simulations 
which show that massive dark matter halos (the very halos in 
which these massive galaxies likely reside) grow by factors 
of > 3 since z = 1 . The success of model Sat2ICL resolves 
this tension by "hiding" the accreted stars in the ICL. Obser- 
vations of the evolution of the ICL at z < 1 will be needed to 
substantiate this picture. 

The success of model Sat2ICL provides us with further in- 
sight into the nature of the ICL. Observationally the ICL ap- 
pears to have colors consistent with the B CG and thus con- 
tains primarily of old stars formed at z > 1 dKrick et al.ll2006t 
iGonzalez etal]l2000t [Zibetti et al.1l2005h . Is this consistent 
with model Sat2ICL, where the ICL is built up by mergers 
with the BCG at z < 1 ? The answer to this question is most 
likely yes, because the subhalos that are disrupting at z < 1 
were accreted onto the host halo at z ~ 1 (the average accre- 
tion epoch, weighted by the fraction of stellar mass brought 
in by the subhalo is 0.93). In other words, the galaxies that 
are contributing stars to the ICL at late times were part of the 
main halo at z ~ 1 and hence could reasonably have had their 
star formation truncated by one or more cluster-specific pro- 
cesses (e.g. ram pressure stripping or harassment) by z ~ 1. 
Thus these galaxies would be adding primarily old red stars 
to the ICL when they disrupt. 

There are several other observables which can provide ad- 
ditional tests of these models. Specifically, the number of 
galaxies, N(M), and the total cluster luminosity, L tot (M), both 
a function of cluster virial mass, provide independent con- 
straints compared to the observations explored herein. How- 
ever, these two observables are much more sensitive to star 
formation since z = 1 (because they include lower mass galax- 
ies), which has been neglected in these models. Hence in the 
present work we have not not included a comparison to these 
observables because such a comparison would require addi- 
tional, less constrained assumptions. 

5.2. Comparison to Related Work 

Semi-analytic models governing the formation and evolu- 
tion of galaxies have proven capable of reproducing both 
strong and mild e volution of massive galaxi es since z = 
1 . The mod els of iKitzbi chler & White! d2007l KW06) and 
iBower et al.l (120061) produce relatively mild evolution in the 
massive end of the GSMF, although KW06 appear to over- 
predict Jhe_jy3undanc^_c^ m gal axies at z = 0. Mean- 
while, [De Lucia & Blaizotl d2007l) and iDe Lucia et all (120061) 
use a semi-analytic model very similar to that used in KW06 
and found that massive galaxies roughly double in stellar mass 
since z = 1 . This doubling in stellar mass does not strongly af- 
fect the evolution in the GSMF presented in KW06 because 
these authors assume 0.25 dex uncertainty in the observed 
stellar mass estimates at z > 0. This assumed uncertainty, 
which is likely an upper bound, increases the model GSMF 
at the massive end for z > 0, and hence masks the stronger in- 



trinsic evolution in the model. None of these models attempt 
to model the ICL, and all are quite sensitive to their treatment 
of the merging of satellite galaxies (the various possible treat- 
ments are not explored in these models) as well as an array 
of model parameters. For these reasons it is difficult to draw 
general conclusions from these models. 

Our approach most closely parallels that of iMonaco et alj 
(2006) who used a semi-analytic model to follow the evolu- 
tion of the GSMF. These authors artificially turned off their 
star-formation prescription in order to follow the dissipation- 
less growth of galaxies since z=l, similar to what we do here. 
When a satellite merges with a central galaxy, they transfer a 
fraction, / sca ttei, of the satellites stars to the ICL. They found 
that /scatter > 0.3 resulted in evolution in the GSMF in agree- 
ment with observations. 

Several models presented in the present work are clo sely 
related to the scheme employed in Mo naco et al.l d2006l) . In 
particular, our models Sat2Cen, Sat2Cen+ICL, and Sat2ICL 
are similar to the ir m odel with / sca ttei = 0.0, 0.5 and 1.0, re- 
spectively 10 . In £ )4.1I we showed that models Sat2Cen+ICL 
and Sat2ICL were indeed in agreement with the observed evo- 
lution of the GSMF while mode l Sa t2C en was not, similar to 
t he c onclusions of M onaco et alJd200 6). However, in £|4.2l and 
£|4.3l we showed that model Sat2Cen+ICL overproduced BCG 
luminosities at z = and underproduced the fraction of com- 
bined BCG and ICL light contained in the ICL, while model 
Sat2ICL successfully reproduced both of these observations. 
We he nce expect that the model presented in IMonaco et alj 
(2006) would reproduce these latter two observations only if 
they used /scatter ~ 1 ■ Comparing models to BCG luminosities 
and ICL fractions vs. cluster virial mass provides unique con- 
straints relative to evolution in the GSMF because these two 
observables directly probe the properties of the most massive 
systems, while the high mass end of the GSMF is sensitive to 
Poisson uncertainty and cosmic variance. 

Each of the models presented herein make predictions for 
the disruption rate of satellite galaxies. IWhite etafl (120071) 
used the redshift-dependent clustering of galaxies from z — 
0.9 to z — 0.5 to constrain the disruption rate of satellite galax- 
ies. Over this time interval, these authors found that at least 
1-2 satellites brighter than > 1.6L* per massive halo were 
disrupted. When focusing on satellites comparable to those 
in White et al., we find on average 1.1 satellites within mas- 
sive halos have disrupted between z = 1 and z = 0.5 (for mod- 
els other than model KeepSat, since in that model satellites 
never disrupt). Agreement between these results provides a 
satisfying cross-check for both approaches. On the one hand, 
we follow directly the evolution of subhalos in simulations 
and hence have useful information regarding their disruption. 
Conversely, White et al. rely primarily on evolution of the ob- 
served clustering of galaxies to constrain the disruption rate, 
and hence make no assumptions regarding the relation be- 
tween subhalos and satellites. When combined, these results 
provide further weight to the idea that, over the mass ranges 
explored herein, disrupted subhalos correspond to disrupted 
satellite galaxies. 

Most recently, Pur cell et al.l ([2007) have modeled the build- 
up of the ICL by tracing the fate of satellites in dark mat- 
ter merger trees generated from the Extended Press-Schechter 
formalism combined with a detailed prescription for satel- 

10 Note that in their model 10% of the stars in the ICL came from the 
tidal stripping of satellite galaxies — a process which we have ignored in the 
present work. 
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lite stripping IZentner et al.l (|2Q05). This dark matter treat- 
ment is coupled with a simple stellar mass-dark matter halo 
correspondence at z = and a simple prescription for star- 
formation. The ICL in their model is built up from dis- 
rupted satellites; with a reasonable threshold for subhalo dis- 
ruption, these authors found trends similar to those reported 
here (since their merger trees are based on analytic prescrip- 
tions, they were able to extend their analysis to halos at con- 
siderably lower mass then explored herein). In particular, they 
found that the fraction of BCG+ICL light contained in the ICL 
in massive h alos is considerable and in agreement with the ob- 
servations of iGonzalez et al.l d2005). The agreement between 
their scheme and ours is encouraging given the substantial dif- 
ferences in the detailed treatment. 

6. SUMMARY 

In this paper we investigated models for the dissipationless 
buildup of massive central galaxies and the intracluster light, 
in the context of merging histories for dark matter halos in 
high resolution simulations of the currently favored ACDM 
model. We used a simple model for associating galaxies with 
dark matter halos and subhalos at z = 1, using the observed 
galaxy stellar mass function and the mass function for ha- 
los and subhalos measured in simulations. The dissipation- 
less evolution of galaxies in this model was tracked with the 
merging history of the dark matter halos extracted from sim- 
ulations. We then confronted this model with data on the evo- 
lution of the galaxy stellar mass function and with the amount 
and fraction of cluster light that is in brightest cluster galaxies 
(BCGs) and in the ICL at z = 0, investigating where the pre- 
dictions from variations in the fate of stars in merging galaxies 
could be distinguished. 

We found that our model accurately reproduces a variety of 
observed properties at z = if disrupted subhalos deposit most 
of their stars into the intracluster light (ICL). Other scenar- 
ios, either those in which most of the stars are deposited onto 
the central BCG or in which stars from disrupted halos are 
left behind as satellite galaxies, are strongly disfavored by the 
data. Such a scenario suggests that, while BCGs do not appear 
to evolve strongly at z < 1 , the ICL surrounding such galax- 
ies is growing substantially over this epoch. This scenario is 
corroborated by high-reso lution dissipationless simula tions of 
galaxy-galaxy mergers (Bovlan-Kolchin & Ma 2007) which 
find that disrupted satellites preferentially build-up the outer 
envelope of massive galaxies. 

Although ideally one should distinguish between light 
bound to galaxies and light that is dynamically bound only 
to the main halo, it is worth noting that our analysis does not 
explicitly distinguish between the outer parts of bright or cD 



galaxies and the ICL. For the purposes of this work "BCG" 
refers to that part of the central galaxy's light that is captured 
in standard survey photometry, operationally defined for most 
of the comparisons herein as light above a surface brightness 
cut of /!,- w 23 mag arcsec" 2 , while "ICL" refers to all cluster 
light that is not contained in galaxies using these measure- 
ments. Further work both on the theoretical and observation 
side is needed to refine this distinction, and care should be 
taken when comparing various studies, as there is a wide vari- 
ation in choices made for these definitions. 

We emphasize that models for the formation and evolution 
of galaxies must be seriously confronted with observations of 
the ICL, in addition to more conventional observations such 
as the GSMF and the two-point correlation function. The ICL 
contains a significant, if still somewhat uncertain, amount of 
stellar mass, and models which ignore this component will 
either place too much stellar mass in resolved galaxies or will 
fail to produce enough stars globally. 

The success of this simple model lends weight to earlier im- 
plications from clustering statistics that the resolution of the 
current generation of iV-body simulations is sufficient to re- 
solve the bulk of subhalos that correspond to observed satel- 
lite galaxies in clusters. Such a confirmation unleashes an ex- 
citing array of possibilities for understanding the connection 
between galaxies and dark matter halos. 
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